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Abstract. 1 A consistent description of the 7t°7Z° invariant mass distribution in the 0(1020) — ► 
/o(980)y — > 7Z°7Z°Y decay and the ^°tj in 0(1020) — > ao(980)y ^ 7Z°rjy is suggested. A search 
for the consequences of the flavor 51/(3) symmetry for the scalar mesons can be based on such an 
analysis. In order to accurately treat the pseudoscalar meson dynamics, which is very important for 
the scalar meson decays, we employ Resonance Chiral Theory. 
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INTRODUCTION 

This paper focuses on the isotriplet ao(980) and isosinglet /o(980) scalar mesons (S). 
These particles, probably together with the controversial o meson, represent the light- 
est members of the scalar meson spectrum. To test this assumption it is important to 
search for consequences of the flavor SI/ (3) symmetry for these particles in nature. 
Among different observables, the invariant mass distributions of %% and %r\ pairs in the 
0(1020) — > Ttity (or nr)y) decays are of considerable interest, because the ao(980) and 
/o(980) are important intermediate states in these decays and thereby show up in the 
spectra. The e + e~ experiments in Novosibirsk [1, 2, 3, 4, 5] and Frascati [6, 7, 8, 9, 10] 
allow one to study the decays. From the point of view of data extraction, the neutral 
final states are preferable since in the decays 0(1020) — ► 7T 7T y, Tt Q r\y the photon may 
be radiated only from the final state. Therefore the above processes are rather suitable 
for study of chiral dynamics. 

The analysis of the data (see Refs. [2, 11] and others) shows that 

• the decays are mediated by the S resonances, /o(980) and ao(980), in the isoscalar 
and isovector channels, respectively; 

• the kaon loop (KL) coupling of S to the vector meson is very important. 

Fig.l shows schematically the processes in question. The KL coupling scheme allows 
one to reproduce a cusp behavior of the amplitude dependence on the invariant en- 
ergy [12]. Many authors relate the dominance of the KL mechanism to a large KK com- 
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FIGURE 1. The e + e~ annihilation to nity (and %T] y). S denotes the intermediate scalar resonance. 



ponent in the ao(980) and /o(980) mesons and to the proximity of the KK threshold to 
the scalar meson mass. We argue that the KL mechanism is a feature of the chiral dy- 
namics and reflects the important role of pseudoscalar mesons in low- and intermediate- 
energy interactions. We also emphasize that the kaon loops in the § — > Sj transitions 
naturally arise in the leading order in Resonance Chiral Theory (R%T), irrespectively of 
the threshold or mass position, and the structure of scalars. In addition, R%T accounts 
for an important feature - momentum dependence of the vertex functions. 

We have recently obtained [13] a complete set of <f?(p 4 ) contributions to various ra- 
diative decays with the scalar mesons from the C{p A ) R%T Lagrangian. The parameters 
were fixed from the average values of the decay widths. In the following we discuss 
the TtyT-based study of the n°n° and K Q r\ pair mass distributions in the 0(1020) — > Sy 
decays and the assumptions used. We fix the parameters and compare predictions with 
the data and other models. 

We conclude that the invariant-mass distribution is a more reliable source of informa- 
tion than the average width from PDG, because the width values may depend on models 
used in the data analysis. 



FORMALISM 

The R%T Lagrangian terms [14] 

L sca lar = C d (S° C %U») + C m (S° Ct X+ ) + C d S S ^ (u^) + C m S S ^ ( X+ ) , (D 
Lvector = ^<WH + ^<^v^ V > (2) 

describe interaction of the scalar singlet S smg , octet S oct , pseudoscalar mesons parame- 
terized by Ufi, vector mesons V^ v and the external electro-magnetic field at order <ff{p A ) 
(for notation we refer to the original paper [14]). The couplings for the scalar singlet c d 
and c m are in general independent of the octet couplings q and c m . The values of Fy and 
Gy are accurately determined from different observables. The flavor SU (3) symmetry is 
explicitly broken by the mass term contained in % + and the relevant couplings have the 
subscript m: c m and c m . We mention below the important aspects of the model [13]. 

• In order to estimate the ciq — > nrj decay width, one has to construct the physical r\ 
state from the pseudoscalar singlet r\\ and the eighth component r/s of the octet. We 
use the two-parameter mixing scheme: r\ = T^cosOs — sin0i with 0$ = —9.2° 
and 0i = -21.2° [15]. 
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FIGURE 2. The invariant mass distributions in <j) — ► n°n°y (left) and — * n°f]y (right). The solid lines 
correspond to the choice of parameters in [13]. The dashed lines are calculated with "classical" values 
for the scalar meson couplings [14] (c m = 42 MeV, c ( [ = 32 MeV, 6 = —35.26"). The dotted lines are 
calculated in the present work. Data: [2] (left) and [5] (right). 



• For the isoscalar state /o(980) we use mixing with the angle 9: fo = S sing cos — 
5g Cf sin0, and for the isovector state ao(980) we take ciq = S% a . This scheme is 
rather flexible and allows us to study the mixing within the lightest nonet. 

• Although the usual large N c behavior of a resonance may not hold for some physical 
scalar mesons (see discussion in Refs. [16, 17]), the scalars S smg and S oct in the 
present model obey this behavior, as these resonances are introduced similarly to 
the vector and axial- vector resonances. Thus we employ the relations c m = c m / y/3 
and Cd = Cd/V3 between the singlet and octet couplings. 

• The effect of the resonance finite width in the invariant mass distributions is impor- 
tant [18]. We choose the propagator of the scalar meson with the mass ms in the 
form 

D S (P 2 ) = \p 2 -m1 + im s r s , tot {p 2 )}-\ (3) 

where the total width T$ tot(p ) depends on the meson momentum squared p 2 . It 
may also be important to use a more advanced form of the propagator including 
both real and imaginary parts of the self-energy, as suggested in Ref. [19]. 

Following approach of Ref. [13] one arrives at the matrix elements for the radiative 
decays. 



DISCUSSION 



Let us summarize advantages of the present approach. 

• There are no contact VyS and yyS couplings at the lowest (tree-level) order. There- 
fore the leading contribution includes one-loop diagrams with intermediate pseu- 
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TABLE 1. The width estimates and comparison with other models 



observable 


estim. 


ri3i* 

1 * *- 1 


[221 * 

1 — — | 


model "a" 
Ref. [23] 


model "b" 
Ref. [23] 


[241** 

L J 


PDG [21] 


r V r " x 10- 4 


0.63 


1.7 


1.4 








0.76 ±0.06 


r 

-V 22 - x 1(T 4 


3.16 


4.4* 


1.4 


4.92 ±0.07 


4.92 ±0.07 




4.40 ±0.21 


V 


4.99 


2.6 


1 


0.26 ±0.06 


0.46 ±0.09 


— 


6.1 ±0.6 


rv^.MeV 

T h ^ K , MeV 


21.11 
54.56 


14.2 
41.8 










T.A. 9+22.7 


TflQ^yy, keV 
r /o ^ rr , keV 


0.16 
0.13 


0.30* 
0.31* 


0.24 
0.24 


0.28 ±0.09* 
0.39±0.13* 


0.28 ±0.09* 
0.39±0.13* 




0.30±0.10 
31 + 08 


Tflo^yp' keV 
r /o ^ rp ,keV 

TaQ-^yfl), keV 

r /f) ^ r<B , keV 


5 
4.1 
4.8 
8.4 


9.1 
9.6* 
8.7 
15.0 


3.4 
3.4 
3.4 
3.4 


3.0±1.0 
4.2± 1.1 
641 ±87 
4.3 ±1.3 


3.0±1.0 
19±5 
641 ±87 
126 ±20 


11 ±4 
3.3±2.0 
31 ± 13 
88 ±17 





* an asterisk marks the input values 

* the kaon loop model estimates are selected from Ref. [22] 

** the model of Ref. [24] accounts also for vector mesons in the loops 

* both pion and kaon loops are included in the current model 



doscalar mesons. This holds for any scalar meson mass irrespectively of their inter- 
nal structure. 

• The sum of the loop contributions is convergent, gauge invariant and universal: one 
can use either the analytical expression of [20], or calculate it numerically. 

One should also stress an influence of the pseudoscalar meson dynamics on the above 
decays, as interaction of the scalar mesons is intimately connected with pseudoscalar 
loops. 

The predictive possibilities of the model are illustrated in Fig. 2. The observables 
shown in the figure are the invariant mass distributions for the neutral pseudoscalar 
meson pairs in the radiative decays. We do not make a numerical fit to the data and 
just illustrate flexibility of the model. The masses of the scalar mesons in the analysis are 
taken mf = 980.0 MeV and m aQ = 984.7 MeV. The parameters are chosen to reproduce 
the form of the distribution in the vicinity of the peak (see the dotted line). Such a 
choice, c m = -38.32 MeV, c d = -9.47 MeV and = -11.62°, also preserves the 
reasonable agreement with the decay widths, see Table 1. The table also shows our 
previous results [13], which were based on analysis of the average width values given 
by PDG [21]. A comparison with predictions of other models [22, 23, 24] is presented 
in Table 1 as well. We also update the predictions for the S — > yV because of the growing 
interest to these decays [13, 22, 23, 24]. One can see that the current approach describes 
the hadronic decays more accurately than the two-photon decays. 

We should notice that the other resonance contributions, which interfere with the 
scalar meson mechanism, are quite complicated, even in case of the neutral particles in 
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the final state (e.g., there are — > p u n u — > JFTry, — > cor/ — > 17 Try, — > p u Tr — > 17 Try 
and other contributions). In this connection we note that a problem of the distribution 
data analysis is discussed in Ref. [25]. It is a challenge to calculate consistently the in- 
terfering contributions in framework of R%T and build a unified analysis of the available 
data. 

The model used here for the mass distributions in the decays can be applied to 
a detailed study of other radiative decays involving the light scalar mesons. Among 
the most interesting ones we should mention the S — > yV decays currently studied 
experimentally in Jiilich [26, 27]. An option to apply the model for the above decays 
and 0(1020) — > TT + TT y data is to include it in a Monte Carlo generator for the Dalitz 
plot or mass distribution analysis. 
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